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Abstract 

We study the synchronization of time-delayed nonidentical maps subject to uni- 
directional (master-slave) coupling. The individual dynamics of the maps have a 
delay n±, and the coupling acts with a delay ri2- We show analytically that, suitably 
tuning the slave map parameters, two distinct synchronization regimes can occur. 
In one regime the lag time between the slave and the master maps is given by the 
delay of the coupling, ri2, while in the other regime is given by the difference be- 
tween the delays, n% — n2- We analyze the effect of the coupling strength on the 
different synchronization regimes in logistic and Henon maps. 
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1 Introduction 



Delay differential equations have received much attention over the years be- 
cause of the significant role of delayed feedback in the dynamics of many 
physical and biological systems [1]. On the one hand, delay-differential sys- 
tems often exhibit multistability -i.e., the coexistence of several attractors- 
and multistability enables such systems to act as memory devices [2,3], an idea 
first suggested by Ikeda and Matsumoto [4]. On the other, the study of delay 
systems is motivated by the fact that these systems exhibit high-dimensional 
chaos and therefore can be used in communication systems based on chaotic 
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synchronization, to securely encrypt information into their chaotic outputs 
[5,6]. 



Several authors have recently shown the existence of two different regimes of 
synchronization in nonidentical coupled time-delayed differential equations, in 
a master-slave configuration [7-11]. When the master and the slave systems 
have both the same amount of delayed intrinsic feedback and the delay time 
r is the same for both systems, the slave system variables synchronize with 
the master system variables with a lag time that is equal to the delay of the 
coupling, t c . In other words, x 2 (t) = Xi(t — T c ), where Xi(t) and X2(t) represent 
the states of the master and slave systems respectively. 

On the other hand, when the master and the slave systems have the same 
amount of delayed intrinsic feedback and the feedback coefficient of the master 
system is equal to the sum of the feedback coefficient of the slave system 
and the coupling coefficient, the systems synchronize with a different lag- 
time. In this case X2(t) = x±(t + r — r c ). Two synchronization regimes are 
therefore possible for appropriate choices of the slave map parameters, as 
shown analytically in specific instances [11,12]. 

In this paper we show that these two synchronization regimes can also occur in 
time-discrete dynamical systems. In Section 2 we consider generic nonidentical 
delayed maps, unidirectionally coupled, and define the two synchronization 
regimes. As an illustration, in Section 3 we analyze logistic and Henon delayed 
maps. Finally, Section 4 presents our conclusions. 



2 Master-slave coupled delayed maps 



We consider a generic master map of the form 

x n +i = af(x n ) + Pf(x n - ni ) + g(x n ), (1) 



where a and (3 are parameters. The slave map is given by 

y n +i = ot s f(y n ) + (3 s f(y n—n\ ) + g{y n ) +vf(x )• (2) 



If the parameters of the slave map are tuned in such a way that a s = a and 
Ps = P — V (case I), the slave map reduces to 

y n +i = af(y n ) + Pf(y n - ni ) + g{y n ) + v[f( x )-f(y n—ni )]■ (3) 
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Full synchronization can be expected for sufficiently large 77 on the synchro- 
nization manifold y n = x n+ni - n2 . On the other hand, if the parameters of the 
slave map are tuned in such a way that a s = a — 77 and f3 s — (5 (case II), the 
slave map reads 

Vn+l = af(y n ) + Pf(y n - ni ) + g(Vn) + V[f( x n-n 2 ) ~ fiVn)], (4) 



and full synchronization may occur on the synchronization manifold y n = 
x n _ n2 . Therefore, depending on the parameters of the slave map, a s and (3 S , 
full synchronization can take place with two different lag times, Aj = n 2 — n\ 
in the case of Eq. (3) and An = n 2 in the case of Eq. (4). Note that, in case 
I, one can have anticipated syncronization for ri\ > n 2 [refs. de anticipated 
synchronization] . 

The actual possibility of observing full synchronization in either case is deter- 
mined by the stability of the synchronized state. Linear stability analysis of 
Eqs. (3) and (4) can be carried out by noticing first that the two equations 
can be written in a unified form as 

n—n\ 71 — 712 )-f(y n-n 3 )], (5) 



with h(y) = otf(y) + g(y)- In Eq. (3) we have = rii, whereas in Eq. (4) 
we have n 3 = 0. The synchronization manifold is y n = x n+n3 _„ 2 . Applying a 
perturbation y n = x n+nz - n2 +S n , replacing in Eq. (5), and taking into account 
Eq. (1) we get, to the first order in the perturbation, 

^n+l h (x n J rnz — n2 S )8 n ~\- (3 f (iC?i+?i3— 7i 2 — ni)^n— ni f]f (^-n— ri2 )^n— n.3 1 (6) 



where primes indicate derivatives. Equation (6) can be formally integrated by 
introducing a linear (N + l)-dimensional map, with N = max{ni,ri3}, for a 
variable r n = (r°, r^, . . . , r%), with r\ = 5 n -k- This equivalent map is given by 

r n+ i = M n r n , (7) 



where the elements of the matrix M n are given by the (time-dependent) coef- 
ficients in Eq. (6) [13]. The solution to Eq. (7) reads 

r n = U n r Q = M n _iM n _ 2 • • • M.MqTo, (8) 



so that the state of full synchronization is linearly stable if all the eigenvalues 
of the evolution matrix U n vanish for n — ■> oo. Whether this condition holds 
or not for a given value of the coupling constant i] can be readily verified by 
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Fig. 1. Schematic representation in the parameter space (a s , (3 S ) of the regions where 
synchronization with lag time Ai and An occur. 

numerical means. Note that all the elements of matrix M n , given by the coef- 
ficients in Eq. (6), involve a delay ri2 which thus acts as a uniform time shift. 
This fact implies that, in the limit n — > oo, the eigenvalues of U n become inde- 
pendent of Tii- Consequently, the value of n<i is irrelevant to the stability of full 
synchronization (cf. [13]). Note carefully that, generally, full synchronization 
for the two cases considered above will be stable on two different ranges of 
the coupling constant 77. If the master map is chaotic, we expect that for suf- 
ficiently small and large r\ full synchronization is unstable and stable in both 
cases, respectively, while for intermediate values only one of the cases admits 
full synchronization. 

We have seen above that full synchronization is possible at two points in the 
(a s ,(3 s ) parameter space, with different lag times in each case. To encounter 
such synchronized states, the slave-map parameters must be exactly tuned 
on one of those synchronization points. Their location is schematically shown 
in Fig. 1. Though when the slave system is slightly detuned with respect to 
the synchronization points full synchronization will not occur, it is expected 
that the slave-map orbit follows approximately the master-map orbit with 
the same lag time, Aj or An- To quantitatively characterize the degree of 
synchronization between the two orbits and the respective lag time, we may 
use the so-called similarity function defined as [11] 

c2 __ ([ X n+A — Un] 2 } 
A [«(^>] 1/2 ' ^ 

where the brakets (•) stand for time averages over asymptotically large times. 
If x n and y n are independent time series with similar mean value and dispersion 
we have <Sa ~ v2 ~ 1.4. If, on the other hand, there is full synchronization 
with lag time A, = 0. The similarity function Sa can be determined, at 
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each point (a s ,f3 s ) parameter space and for each lag time A. At (a, (3 — rj) we 
should have = 0, while at (a—rj, (3) we should have Sa h = 0. It is expected, 
moreover, that in a region around each synchronization point the similarity 
function attains a minimum, as a function of the lag time, for A = Aj and 
A = An, respectively. These regions are qualitatively illustrated in Fig. 1. In 
the remaining of the parameter space, as far as the slave-map orbits do not 
diverge, the similarity function can attain a minimum for any other value of A 
-without reaching, however, = 0. Note that the boundaries of such regions 
will in general depend on the coupling constant rj. In the following sections, 
we study these aspects of synchronization in Eqs. (1) and (2) for logistic and 
Henon delay maps in their chaotic regime. 



3 Synchronization of delayed logistic and Henon maps 

3.1 Logistic maps 

As a first illustration of the synchronization properties of Eqs. (1) and (2) in 
cases I and II, we consider the choice f(x) = x(l — x) and g(x) = 0. The 
master system becomes a delayed logistic map 

x n _|_x <yx n (l Xn) -\- f3x n — ni {\ x n — ni ), (10) 

whose orbits are bounded to the interval (0, 1) for < a, (3 and a + (3 < 4. 
In different regions of the parameter space (a, (3) and depending on the delay 
rii, this system displays periodic, quasiperiodic and chaotic evolution. The 
corresponding slave map is given by 

y n+ i = a s y n (l - y n ) + f3 s y n - ni {l - y n - ni ) + r]x n - n2 {l - x n _„ 2 ). (11) 

Its orbits are nondivergent for < a s , /3 s ,r) and a s + (3 S + tj < 4. 

Figure 2 displays the synchronization regions in the parameter space (a s , f3 s ) 
and how they vary as the coupling coefficient 7] increases. The master map 
parameters are a — 1.8, (3 — 2.1 and the delay times are ni = 2, n 2 = 3. 
For each pair (a s , (3 S ) we have determined as a function of the lag-time 
A, and detected the value of A for which the similarity function attains its 
minimum min(S'A). The left column of Fig. 2 displays this minimum for three 
values of rj. Light tones represent low values of min(S'A), i.e. high master-slave 
correlation, while darker tones correspond to poor correlation [min(S'A) ~ 1]. 
In the black upper-right region the slave-map orbits diverge. 
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The right column of Fig. 2 displays the lag-time A at which the similarity 
function attains its minimum. The region where Sa is minimal with lag-time 
Ai = n 2 — rii = 3 is represented by the darker gray tone, while the region with 
lag-time A n = n 2 = 3 is represented by the lighter gray tone. White represents 
the parameter region where the minimum value of the similarity function 
occurs for a lag-time which is different from Ai or An- Black represents the 
parameter region where the trajectory of the slave map diverges. 

For low coupling intensity, rj = 0.2, both synchronization regimes are unstable 
[Figs. 2 (a, b)]. The synchronization regions are not well defined and have 
fuzzy boundaries. While the minimum of the similarity function at point I 
occurs at the expected lag-time Ai, the minimum of the similarity function at 
point II occurs at a different lag-time. Notice, in fact, that in Fig. 2 (b) point 
II is in the white region that represents a lag-time different from Aj or A n . 
The region corresponding to each regime is disconnected and quite complex in 
shape, with parts in distant zones of the parameter space. Note, for instance, 
the light-gray zones near (3 S = where master-slave correlation is however 
rather poor. 

As the coupling intensity grows, zones I and II become more uniform and 
increase in total extension. For rj = 0.8 only regime I is stable [Figs. 2 (c, d)]. 
In this case, min(SA) at points I and II is equal to and 0.2 respectively. 
For large enough r], both regimes are stable [Figs. 2 (e, f); r] = 1.2] and, as 
expected, min(SA) equals zero at points I and II. 

Figure 3 illustrates the master-slave correlation at different points of parameter 
space for r\ = 1.0 (all other parameters are as in Fig. 2). In this case, type 
I synchronization is stable [Fig. 3 (a)] but type II is not [Fig. 3 (b)], and 
it is worth mentioning that the minimum value of the similarity function, 
Sa = 0.057, does not occur at point II (a s = 0.8, f3 s = 2.1) but at a point 
close to it (a s = 0.75, f3 s = 2.1). Figure 3 (c) displays the correlation plot at 
a point where the lag-time at which Sa attains its minimum is A = —15. 



3.2 Henon maps 

As a second example, we study now a master-slave configuration where each 
element is a two-dimensional delay map, namely, a Henon-like map. The 
evolution of the master coordinates x„ = (u n , v n ) is given by the functions 
f (x) = (— -u 2 , 0) and g(x) = (1 + v , bu), so that the master system is 

u n+l = l-au 2 n - pu 2 n _ m + v n , ^ 

V n+1 = bu n , 
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Fig. 2. Synchronization regions in the (a s ,(3 s ) parameter space in the case of the 
logistic map, for increasing coupling: (a, b) rj = 0.2; (c, d) r] = 0.8; (e, f) rj = 1.2. 
The master map parameters are a = 1.8, = 2.1, and the delay times are ri\ = 2 
and n-2 = 3. The right column displays the minimum of the similarity function. 
Light tones represent low values of min(S'A) (good master-slave correlation) and 
vice versa. Black represents the region where the slave-map trajectories diverge. 
The left column displays the lag-time where the minimum value of Sa occurs. In 
the light-gray region, the lag-time is An = n2 = 3, while in the dark-gray region it 
is Ai = ri2 — n\ = 1. In the white region the lag-time is different from Ai or An- 
The small circles stand at the synchronization points I and II, (a s = a, (3 S = (3 — r]) 
and (a s = a — r/), S = (3), respectively. 

cf. Eq. (1). In the following we choose b = 0.3. The slave system, with coordi- 
nates y n = (w n , z n ), is governed by the equations 



w n+1 = 1 - a s w 2 n - /3 s w 2 n _ ni + z n - 77<_ n2 , 
z n+ i = bw n , 



(13) 
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Fig. 3. Correlation plots for r] = 1.0, n\ = 2, ri2 = 3 and (a) a s = 1.8, (3 S = 1.1, 
min(SA) = 0; (b) a s = 0.75, 0, = 2.1, min(S , A ) = 0.057; (c) a s = 1.25, (5 S = 1.525, 
min(S'A) = 0.25. The master-map parameters are as in Fig. 2. 



so that coupling acts on the first coordinate only. The synchronization mani- 
fold is given by 



w n = u 



n+ri3—n2 i 



(14) 



z n = v. 



n+n3—n2 ) 



where, as before, n 3 = n x in case I and n 3 = in case II. 

Next we study in which regions of parameters the different synchronization 
regimes occur. We take parameters for the master map a, ft, n\ such that its 
dynamics is chaotic. 

Figure 4 displays the minimum of the similarity function and the lag-time for 
which the minimum occurs, in the parameter space (a s ,/3 s ). The results are 
similar to those found with the logistic map. For weak coupling the synchro- 
nization regions are not well defined, but as the coupling increases their size 
grows and the boundary between them becomes well defined. For large rj both 
synchronization regimes are stable. Thus, it can be though that a small vari- 
ation of the slave map parameters a s or (5 S near the boundary region might 
induce a transition from synchronization with lag-time Ai to synchronization 
with lag-time An or viceversa. However, near the boundary region we find 
min(S'A) ~ 0.5, which indicates bad synchronization. Therefore, while the lag- 
time at which the minimum value of Sa occurs changes abruptly (from Ai to 
An), there is no sharp transition between one regime of synchronization to 
the other. If the slave map parameters are gradually modified from point I to 
point II, synchronization with lag-time Ai is gradually lost, and as we enter 
the region II, synchronization with lag-time An is gradually established. 
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Fig. 4. Synchronization regions for delayed Henon maps. The parameters of the 
master map are a = 0.7, [3 = 0.7, and delay times are n\ = 6 and 722 = 3. The 
left column displays the minimum of the similarity function and the right column 
displays the lag-time where min(S'A) occurs. In the light-gray region, the lag-time 
is An = ri2 = 3, while in the dark-gray region it is Ai = m — n\ = —3. In the white 
region of (b) the lag-time is different from Aj or An- The small circles indicate the 
points (a s = a — 77, (3 S = (3), and (a s = a, j3 s = (3 — 77. (a, b) 77 = 0.2; (c,d) r/ = 0.4; 
(e, f) r] = 0.6. 

4 Conclusion 



We have studied two regimes of synchronization of delayed nonidentical maps. 
We have shown analytically that, by suitably tuning the slave map parameters, 
two distinct synchronization regimes can occur. In one regime the lag time 
between the slave and the master maps is given by the delay of the coupling, 
77, 2 , while in the other regime is given by the difference between the delays, 
n 1 — n 2 . We have also shown that these two regimes are actually two particular 
cases of synchronization with a slave map that is identical to the master map 
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but that has two delayed feedback terms. 

The two synchronization regimes has been exemplified by considering delayed 
logistic and Henon maps. In both cases, the synchronization regimes are si- 
multaneously stable only for large values of the coupling rj, and therefore, they 
occur at parameters of the slave map, (a s ,(3 s ), which are far away from each 
other. In other words, our results show that in the case of delayed logistic and 
Henon maps, a small variation of a parameter of the slave map can not induce 
a transition from regime I to regime II or viceversa, since they occur in distant 
regions of the parameter space. On the contrary, in the case of semiconductor 
lasers with optical feedback, it has been shown numerically [8,10] that close 
to the lasing threshold, by carefully tuning a parameter of the slave laser one 
can induce a transition from one regime of synchronization to the other. It 
will be interesting to study a delayed map which shows this type of transition, 
and allows an analytical investigation of the phenomenon. 
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